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Context. It has recently been shown that the terrestrial planets and asteroid belt can be reproduced if the giant planets underwent an inward- 
I I ■ then-outward migration (the "Grand Tack"; Walsh et al 201 1). Inward migration occurs when lupiter opens a gap and type II migrates inward. 
P^ \ The planets "tack" and migrate outward when Saturn reaches the gap-opening mass and is caught in the 3:2 resonance with lupiter. 
\J^ Aims. The aim is to test the viability of the Grand Tack model and to study the dynamical evolution of lupiter and Saturn during their growth 
(-H from 10 Me cores. 

Q,. Methods. We have performed numerical simulations using a grid-based hydrodynamical code. Most of our simulations assume an isothermal 
I equation of state for the disk but a subset use a fully-radiative version of the code. 

O Results. For an isothermal disk the two phase migration of Jupiter and Saturn is very robust and independent of the mass-growth history of 
-j-J , these planets provided the disk is cool enough. For a radiative disk the we find some outcomes with two phase migrations and others with more 
^ ■ complicated behavior. We construct a simple, 1-D model of an evolving viscous disk to calculate the evolution of the disk's radiative properties: 
I 1 1 the disk transitions from radiative to isothermal from its outermost regions inward in time. 

Conclusions. We show that a two-phase migration is a natural outcome at late times even under the limiting assumption that isothermal 
^~~^ conditions are required. Thus, our simulations provide strong support for the Grand Tack scenario. 
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I I [ 1. Introduction 3. Once fully-formed, Saturn migrated faster (Masset & 

O' ^, ,. ^ , . . ... .- , o , Papaloizou 2003), caught up to Jupiter, and was trapped in 

. The current paradigm of the origin and evolution ol the Solar , _ ,„. n xr i ^r,r,ox t ■ i i ■ 

T— I ^ ,.,.-„ , ,. . 3:2 resonance (Pierens & Nelson 2008). Interestingly, this 

, System s giant planets follows several distinct stages: i • i^ j i i r , ■ , 

'""I . -' ° ^ ° result IS found to be a very robust outcome of the simula- 

^ I 1. The cores of Jupiter and Saturn form by accretion of plan- tions, independent on the earlier evolution of Saturn's core. 

etesimals (e.g., Kokubo & Ida 1998; Levison et al. 2010). For instance, Pierens & Nelson (2008) investigated the sce- 

The timescales of accretion are poorly constrained because nario in which Saturn's core is initially trapped at the edge 

j^ . as they grow the cores migrate due to both the back-reaction of Jupiter's gap and grows through gas accretion from the 

■ ■ ' from planetesimal scattering (Fernandez & Ip 1984; Kirsh disk. In that case, they demonstrated that although Saturn 

et al. 2009) and type I (tidal) interactions with the gaseous is temporarily locked in the 2:1 resonance with Jupiter, it 

protoplanetary disk (Goldreich & Tremaine 1980; Ward becomes ultimately trapped in the 3:2 resonance. 

1986; Paardekooper et al. 2010). 4. Once Jupiter and Saturn are trapped in 3:2 resonance, the 

2. Jupiter and Saturn's cores slowly accrete gas and each un- gaps carved by the two planets in the Solar Nebula overlap. 

dergo a phase of rapid gas accretion (e.g., Mizuno 1980, Saturn's gap is not as deep as Jupiter's (due to its smaller 

Pollack et al. 1996). The rapid phase of accretion is trig- mass), and this causes Jupiter and Saturn to migrate out- 

gered when the mass in each planet's gaseous envelope is ward while remaining in 3:2 resonance, provided that both 

comparable to the core mass. Runaway gas accretion lasts the disk thickness and the disk viscosity are small enough 

for roughly the local Kelvin-Helmholtz time and ceases (Masset & Snellgrove 2001; Morbidelli & Crida 2007). 

when the planet opens an annular gap in the disk and tran- Outward migration is stopped when the disk dissipates or, 

sitions to type II migration (Lin & Papaloizou 1986; Ward if the disk is flared, at a critical distance where the disk is 

1997). Because of its larger mass and smaller orbital radius, 

Jupiter is thought to have undergone runaway gas accretion 

before Saturn. 
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too thick and the structure of the two planets' common gap 
is compromised (Crida et al 2009). [J 
5. After the dissipation of the gas disk, planetesimal-driven 
migration causes a large-scale spreading of the planets' or- 
bits because Jupiter is the only planet for which the ejec- 
tion of small bodies is more probable than inward scat- 
tering (Fernandez & Ip 1984; Hahn & Malhotra 1999). In 
the "Nice model", Jupiter and Saturn are assumed to have 
formed interior to their mutual 2:1 resonance and, when 
they cross it, an instability is triggered that causes the Late 
Heavy Bombardment (Gomes et al. 2005). Recent work 
has shown that the Nice model is still valid if more real- 
istic initial conditions are used, with Jupiter and Saturn in 
3:2 resonance and Uranus and Neptune also trapped in a 
resonant chain (Morbidelli et al. 2007; Batygin & Brown 
2010). The Nice model can reproduce the giant planets' fi- 
nal orbits (Tsiganis et al. 2005), the orbital distribution of 
Jupiter's Trojan asteroids (Morbidelli et al. 2005), and sev- 
eral other characteristics of the Solar System's small body 
populations. 

Although the detailed orbital evolution of the giant planets 
is not known, these steps explain their origin in broad strokes. 
By assembling steps 2-4, Walsh et al. (201 1) recently proposed 
a new model to explain the origin of the inner Solar System 
called the "Grand Tack". In this model, Jupiter formed at ~ 2-5 
AU, migrated inward then "tacked" (i.e., changed the direction 
of its migration) at an orbital distance of ~ 1 .5 AU when Saturn 
caught up and was trapped in 3:2 resonance and migrated back 
out past 5 AU. Jupiter's tack at 1 .5 AU truncates the inner disk 
of planetary embryos and planetesimals from which the terres- 
trial planets formed at about 1 AU. This type of narrow trun- 
cated disk represents the only initial conditions known to satis- 
factorily reproduce the terrestrial planets, in particular the small 
mass of Mars compared with Earth (Wetherill 1978; Hansen 
2009; Raymond et al. 2009; Walsh et al. 2011). An additional 
success of the Grand Tack model is that the asteroid belt is nat- 
urally repopulated from two distinct populations corresponding 
to the C- and S- type asteroids. At the end of the Grand Tack, 



' We note that this step is somewhat uncertain because both the mi- 
gration and accretion rates of giant planet cores should a roughly linear 
dependence on the planet mass. Thus, we naively expect that Saturn's 
gas accretion should mimic Jupiter's as it migrates inward, meaning 
that Saturn should be roughly IMj when it catches up to Jupiter, pre- 
cluding outward migration, which requires a Saturn/Jupiter mass ratio 
of roughly 1/2 or smaller (Masset & Snellgrove 2001). The solution 
to this problem is not clear: it may involve a change in the disk opac- 
ity to allow Saturn's rapid type III migration to last for longer than 
Jupiter's. Such rapid migration has been invoked to explain the 3:2 
resonant exoplanet system HD 45364 (Rein et al. 2010). Of course it 
is reasonable to expect that this mechanism is probably not universal, 
since it depends on details such as the timing of core formation and 
the disk properties. Thus, in many exoplanet systems "Saturn" would 
have reached IMj and the two planets would not tacked and migrated 
outward. The architecture of such systems naturally would not resem- 
ble the Solar System. Understanding the statistical distribution of gi- 
ant exoplanetary systems can therefore place constraints on their early 
evolution and the frequency of "grand tacks". 



the giant planets' orbits represent the initial conditions for the 
Nice model (Raymond et al., in prep.). 

The goal of this paper is to test the viability of the evolution 
of Jupiter and Saturn in the Grand Tack model (Walsh et al 
2011). To accomplish this we use the GENESIS hydrocode to 
simulate the growth and migration of Jupiter and Saturn from 
10 Me cores. With respect to previous simulations (Masset & 
Snellgrove 2001; MorbidelU & Crida 2007; Pierens & Nelson 
2008), we consider a self-consistent scenario in which the cores 
of Jupiter and Saturn slowly grow to full-fledged gas giants by 
accreting gas from the disk. We find that a two phase migration 
of Jupiter and Saturn is a very robust outcome in isothermal 
disks, but occurs in only one of two simulations in radiative 
disks. We place our simulations in the context of the evolving 
Solar Nebula using a 1-D diffusion algorithm that differentiates 
between radiative and isothermal behavior. Our results strongly 
favor a two-phase migration of Jupiter and Saturn, and support 
the Grand Tack. 

The paper is organized as follows. In Sect. 2, we describe 
the hydrodynamical model. In Sect. 3, we present the results 
of isothermal simulations. In Sect. 4 we present results of ra- 
diative simulations. In Sect. 5 we construct a 1-D model of the 
Solar Nebula to show when the disk should be isothermal or 
radiative. Finally, we discuss our results and draw conclusions 
in Sect. 6. 



2. The hydrodynamical model 

2.1. Numerical method 

In this paper, we adopt a 2D disk model for which all the 
physical quantities are vertically averaged. We work in a non- 
rotating frame, and adopt cylindrical polar coordinates (R, (p) 
with the origin located at the position of the central star. 
Indirect terms resulting from the fact that this frame is non- 
inertial are incorporated in the equations governing the disk 
evolution (e.g. Nelson et al. 2000). These are solved using 
the GENESIS hydrocode for which a full description can be 
found for example in De Val-Borro et al. (2006). The evolu- 
tion of each planetary orbit is computed using a fifth-order 
Runge-Kutta integrator (Press et al. 1992) and by calculating 
the torques exerted by the disk on each planet. The disk ma- 
terial located inside the Hill sphere of each planet is excluded 
when computing the disk torques (but see §below). We also 
employ a softening parameter b = 0.6// , where H is the disk 
scale height, when calculating the planet potentials. 

In the simulations presented here, we use Nr = 608 ra- 
dial grid cells uniformly distributed between /?,„ = 0.25 and 
Rom = 7 and A^^ = 700 azimuthal grid cells uniformly dis- 
tributed between 0„„„ - and (pi„ax = 27r. Wave-killing zones 
are employed for R < 0.5 and R > 6.5 in order to avoid wave 
reflections at the disk edges (de Val-Borro et al. 2006). 

For most of the simulations, we adopt a locally isothermal 
equation of state with a fixed temperature profile given by T = 
To(R/Roy^ where /3 - I and where To is the temperature at 
/?() = !. This corresponds to a disk with constant aspect ratio h 
for which we consider values of /; - 0.03,0.04, 0.05. The initial 
surface density profile is chosen to be 2(/J) - 2o(^/^o) '^ with 
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(J - 1 12, 3/2 and Eo = 4 x lO"**. In our units, this corresponds 
to a disk containing ~ 0.04 Mq within 40 AU. 

The adopted computational units are such that the mass of 
the central star M, = 1 corresponds to one Solar mass, the 
gravitational constant is G = 1 and the radius R - \ in the 
computational domain corresponds to 5 AU. Thus, in the rest 
of the paper time is measured in units of orbital periods siR - I 
{~ 10 years for 5 AU). However, we note that these simulations 
can be scaled to different disk parameters. For example, if we 
assume R - \ coiTesponds to 1 AU, then the coiTesponding disk 
mass is roughly 10 times larger (for our fiducial iJ"'''^ surface 
density profile). To run simulations that better reproduce the 
conditions at ~ 1 AU would require integration times that are 
ten times longer with disks that are one tenth the density, and 
are not computationally feasible given our current resources. 
Thus, our simulations are intended to demonstrate the relevant 
mechanisms in a similar setting that is, admittedly, somewhat 
more distant. However, we did perform one simulation in a very 
low-mass disk in which R - I corresponded to 1 AU. That 
simulation, presented briefly in Sect. 5, serves to validate our 
results. 

We have also performed a few additional radiative runs 
where the thermal energy is solved. Source terms correspond- 
ing to viscous heating and local radiative cooling from the 
disk surfaces are implemented and handled similarly to Kley & 
Crida (2008), except that we use the Rosseland mean opacity 
given by Bell & Lin (1994). Heat diffusion in the disk midplane 
is not taken into account in the simulations presented here. 

Viscous stresses probably arising from MHD turbulence are 
modelled using the standard 'alpha' prescription for the disk 
viscosity v = aCsH (Shakura & Sunyaev 1973), where c., is the 
isothermal sound speed and H is the disk scale height. In the 
simulations presented here, we use values of a = 2 x 10"^ and 
a = 4x 10 I 

2.2. Initial conditions 

In our simulations, the masses of Jupiter nij and Saturn nis are 
initiated with values of mjj - nigj - 10 M®; and the cores 
of Jupiter and Saturn are placed on circular orbits at oy = 2 
and as - 2.65, just exterior to their mutual 3:2 mean motion 
resonance. 

In an isothermal disk, the type I migration timescale of a 
planet with mass nip, semimajor axis Oj, and on a circular orbit 
with angular frequency Qp can be estimated by (Paardekooper 
etal.2010): 



= (1.6+^-H0.7<r)-'^-%/z^O,-'. 
mp I.(ap)af, 



(1) 



Because of Eq.[Tl we expect Jupiter and Saturn's cores, embed- 
ded in a disk model with cr < 3/2 , to undergo convergent mi- 
gration and become eventually trapped in the 3 : 2 resonance. 
In contrast, a disk model with cr > 3/2 should lead to divergent 
migration. 

Once the cores have evolved for ~ 500 orbits of the in- 
nermost embryo, we allow Jupiter's core to accrete gas from 
the disk. For each timestep At, accretion is modeled by re- 
ducing the surface density in the grid cells located within a 



distance Race of the planet by a factor 1 - fjAt. Following 
Paardekooper & Mellema (2008), we set fj = 5/3 in our sim- 
ulations. Furthermore, we choose Race = 0.1 Rh,j where Rfjj 
is the Hill radius of the planet {Rhj - aj(mj/3MQ)^'^); this 
value is small enough to ensure that the accretion procedure is 
independent of our choice of fj (Tanigawa & Watanabe 2002). 

Accretion onto Saturn's core is handled in the same way but 
we use different values for the accretion parameter fs . In this 
work, this process is switched on when Jupiter has grown to a 
fraction xj - nij/Mj of its final mass Mj = 318 M®. We varied 
the value for xj in such a way that different starting times for 
accretion onto Saturn's core are considered. For instance, we 
used xj - 0.5 in most of the simulations presented here, mean- 
ing that accretion onto Saturn's core is triggered once Jupiter 
has reached half of its final mass. We also tested values for xj 
of (concurrent accretion of both cores) and 1 (isolated growth 
of each core in succession). 

The parameters for all of our locally isothermal simulations 
are shown in Table [T] Our fiducial simulation (model II) had 
h = 0.04, cr = 0.5, Q- = 2 X 10-^ fj = fs = 5/3, and cr = 1/2. 
We performed additional runs varying the disk aspect ratio, vis- 
cosity, and surface density profile, varying each parameter or- 
thogonally to our fiducial case and keeping the other parame- 
ters fixed. 



3. Isothermal simulations 

3.1. Evolution of our fiducial case 

Fig. [T] shows the evolution of our fiducial model II in which 
Saturn starts to accrete gas once Jupiter has grown to half of its 
final mass. The time evolution of the planet masses is shown in 
the first panel. It is worth noting that the gas accretion rate onto 
Saturn's core is actually higher than Jupiter's. This is simply 
because the disk mass included in the feeding zone of a planet 
with semimajor axis ap increases as ~ flp '^. The second panel 
depicts the evolution of the semimajor axes. In this run, the 
two cores migrate convergently at early times and are captured 
in the 3:2 resonance. The resonance is maintained until t ~ 500 
orbits when Jupiter starts to accrete gas from the disk. Then, 
the two cores move away from each other because Jupiter's 
migration rate increases as ~ mj before it opens a gap. This has 
the consequence of disrupting the 3:2 resonance. At later times, 
the migration once again becomes convergent as Jupiter clears 
a gap and transitions to slower, type II migration, while Saturn 
starts accreting gas and its type I migration accelerates. The 
fourth panel of Fig. [T] shows the evolution of the period ratio 
(as /aj)^^^. The slowing down of Jupiter's migration due to the 
onset of non-linear effects (once nij < Q.3Mj) is illustrated by 
the presence of a local maximum in the period ratio at f ~ 1200 
orbits, while the onset of gas accretion onto Saturn's core is 
responsible for the sudden drop in period ratio at f ~ 2000 
orbits. 

Convergent migration causes the 3:2 resonance to be recov- 
ered at f ~ 2300 orbits and the resonance remains stable for the 
duration of the simulation. This is illustrated by the upper panel 
of Fig. |3] which shows the time evolution of the two resonant 
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Table 1. Parameters used in the simulations. 
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Fig. 1. Upper left (first) panel: the time evolution of the masses of Jupiter (solid line) and Saturn (dashed line) for run II . Upper 
right (second) panel: time evolution of the semimajor axes. Third panel: time evolution of the eccentricities. Fourth panel: time 
evolution of the period ratio (aslajY'^. 



angles associated with the 3:2 resonance: 



tp - 3/ls - 2/ly - uj and 



3As-2Aj-us, (2) 



where Aj (As) and wy {<jl>s) are the mean longitude and longi- 
tude of pericentre of Jupiter (and Saturn). Capture in the 3:2 
resonance causes the eccentricities of both planets to grow to 
ej ~ 0.03 and es ~ 0.1, as seen in the third panel in Fig.[T] In 
agreement with previous studies (Masset & Snellgrove 2001, 
Morbidelli & Crida 2007, Pierens & Nelson 2008), the long- 
term outcome for model II after capture in the 3:2 resonance 
is outward migration of the Jupiter-Saturn system with the two 
planets maintaining the 3:2 resonance and sharing a common 



gap. Here, the migration reversal occurs at f ~ 2500 orbits, 
when nij ~ 0.6 Mj and ;«s ~ 0.3 Mj. These values are in 
reasonable agreement with Masset & Snellgrove (2001) who 
estimated a critical mass ratio of mjmj < 0.62 for the posi- 
tive torque exerted by the inner disk on Jupiter to be larger than 
the negative torque exerted by the outer disk on Saturn. Fig.|2] 
shows a snapshot of the disk at a point in time where Jupiter 
and Saturn are fully-formed, locked in 3:2 MMR and migrate 
outward. 

For this model, we observe a trend for the amplitude of the 
resonant angles to slightly increase with time to such an extent 



Pierens & Raymond: Two-phase migration of Jupiter and Saturn 



> 




1.000 



Fig. 2. Snapshot of the perturbed disk surface density for model 
II at a point in time where Jupiter and Saturn are fully formed, 
locked in a 3:2 MMR and migrate outward. 



that (p switches from libration to circulation at f ~ 1.3 x 10^ 
orbits (see upper panel of Fig. [3]l. This causes not only the out- 
ward migration rate to subsequently slow down (second panel) 
but also the eccentricities to slightly decrease to values such 
that ej ~ 0.02 and es ~ 0.04 at the end of the simulation (third 
panel). It is interesting to note that some of the known exo- 
planet systems appear to exhibit libration of only one resonant 
angle. In the case of HD128311 (Vogt et al 2005), this behavior 
might be explained by a scattering event (Sandor & Kley 2006) 
or by effects due to turbulence (Rein & Papaloizou 2009). 

In simulation II, the change from libration to circulation 
of one resonant angle occurs when Saturn is at ~ 3.5 numeri- 
cal units, roughly half the distance of the outer disk edge. To 
test whether this could have been a numerical effect caused by 
unresolved high-order density waves (perhaps associated with 
the outer 3:1 resonance with Saturn) we re-ran the latter evo- 
lution of simulation II but with an outer disk boundary at 5 
numerical units rather than 7. In the test simulation, the change 
from libration to circulation again occurred when Saturn was 
at ~ 3.5 units, showing that the outer disk was not the cause of 
the change. 

It appears that one possibility that caused the shift from 
libration to circulation of the resonant angle was the corota- 
tion torque exerted on Saturn. Examination of the torques ex- 
erted on Saturn indeed reveals a tendency for the torque ex- 
erted by the disk material located in between the two planets to 
increase. This suggests that, as the orbital separation between 
Jupiter and Saturn increases during their outward migration in 
3:2 resonance, gas flowing across Saturn's orbit can feed the 
common gap and consequently exert a positive horseshoe drag 
on Saturn. Such a process tends to further push Saturn outward 
(Zhang & Zhou 2010), which may subsequently lead to the 
disruption of the apsidal corotation. This effect may have im- 
portant consequences for the ability of two planets to undergo 



long-range outward migration via the Masset & Snellgrove 
(2001) mechanism. 

3.2. Effect of simulation parameters 

In the following, we investigate how the evolution of the sys- 
tem depends on the simulations parameters. We test the effect 
of varying the start time of Saturn's growth, disk aspect ratio, 
viscosity and surface density profile. Also not shown here, we 
have also tested the effect of changing the accretion parameter 
fs of Saturn and have performed one simulation with fs = 5/6 
(model 12). In that case, the evolution of the system was found 
to be very similar to that obtained in model II. 

3.2.1 . Dependence on the start time of Saturn's 
growth 

In order to examine how the evolution depends on the mass- 
growth history of Jupiter and Saturn, we performed two addi- 
tional simulations varying the time when Saturn's core starts 
to accrete gas from the disc (Table 1). In simulation 13, accre- 
tion onto the cores of Jupiter and Saturn are switched on at 
the same time while in simulation 14 Saturn's accretion started 
once Jupiter was fully formed. Fig. |4] shows the disk surface 
density profiles, the positions of the planets and the locations 
of the 2:1 and 3:2 resonances with Jupiter for simulations II, 
13, and 14 just before accretion onto Saturn's core is switched 
on As Jupiter grows, we see a clear tendency for Saturn's core 
to follow the edge of Jupiter's gap where its inward migration 
(caused by its differential Linblad torque) is balanced by the 
corotation torque (Masset et al. 2006). Of particular importance 
is the location of the gap edge with respect to Jupiter In runs II 
and 14 the edge of Jupiter's gap is located just outside the 3:2 
resonance whereas for run 13 it lies beyond the 2:1 resonance. 

Fig. |5] shows the evolution of simulations with different 
start times for Saturn's accretion {xj = 0,0.5, 1). It is inter- 
esting to note that Saturn's accretion acts to accelerate Jupiter's 
accretion (see upper panel). Jupiter grows to its full mass in 
~ 2900 and ~ 1900 orbits in models II and 13, respectively, but 
takes ~ 7000 orbits to acquire its final mass in model 14. This 
occurs because, as Saturn grows and begins to form a gap, the 
disk surface density near Jupiter increases, thereby enhancing 
Jupiter's accretion rate. 

The second panel in Fig. |5] shows the orbital evolution of 
Jupiter and Saturn. The three simulations behave similarly be- 
fore Jupiter starts to accrete gas from the disk at f ~ 500 orbits. 
All cases undergo convergent migration of both cores followed 
by capture in the 3:2 resonance. At later times, however, the 
evolution of the three simulations diverge. For model 13, the 
fact that gas accretion onto Saturn's core proceeds more rapidly 
compared with Jupiter leads to an even faster convergent migra- 
tion compared with model II, driving the planets even deeper 
into 3:2 resonance after resonant locking. This is illustrated in 
the second panel of Fig.[3]which shows, for this simulation, the 
evolution of the resonant angles associated with the 3:2 res- 
onance. In contrast with run II, periods of circulation of the 
resonant angles are not observed. Instead, gas accretion onto 
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Fig. 3. For models II, 13 and 14 the time evolution of the reso- 
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the cores make the libration amplitudes of the resonant angles 
decrease for t < 2000 orbits. 

As before, once Jupiter and Saturn have grown to nij ~ 0.5 
Mj and nis ~ 0.3 Mj, both planets migrate outward in concert, 
this time at a slightly faster rate than model II. This is because, 
as discussed above, in model 13 the faster convergent migration 
at ? < 2000 orbits due to Saturn's growth has the consequence 
of locking the planets more deeply in resonance than in model 

11 ( i.e., the libration width of the resonant angles is smaller for 

12 than II; see Fig.[3]l. At later times however, the amplitude 
of the resonant angles increases and outward migration slows. 
Thus, the final outcome for run 12 is very similar to that of 
model II in terms of both the late-time outward migration rate 
and the planets' eccentricities. 
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Fig. 5. Upper (first) panel: the evolution of the planet masses 
for models II (black), 13 (red) and 14 (blue). Second panel: evo- 
lution of the semimajor axes. Third panel: evolution of Jupiter's 
eccentricity. Fourth panel: evolution of Saturn's eccentricity. 
Fifth panel: evolution of the period ratio. 
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Fig. 6. Evolution of the torques exerted on Jupiter (left panel) 
and Saturn (right panel) for models II, 13 and 14 at times where 
a quasi-stationary state is reached. 




Fig. 7. The disk surface density profile for models II (black), 
13 (red) and 14 (blue) at f ~ 1.5 x 10"* orbits. 
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For simulation 14 - in which gas accretion onto Saturn's 
core started only after Jupiter reached its final mass - the evo- 
lution differed significantly from runs II and 13 (Fig. |5]l. Up 
to f ~ 2000 orbits, which corresponds to nij ~ 0.5 Mj, the 
evolution of the system is similar to model II. After this point, 
however, the evolution of run 14 diverges from run II. Jupiter 
type II migrates inward whereas Saturn's core migrates slightly 
outward (second panel in Fig.|5]l. This is because Saturn's core 
follows the edge of Jupiter's gap where the positive corota- 
tion torque balances the negative differential Lindblad torque 
(Masset et al. 2006). As Jupiter grows, its gap slowly widens 
causing the orbital separation between the two planets to in- 
crease beyond the 2 : 1 resonance (Fig. nil. Once Jupiter reaches 
its final mass, at f ~ 7500 orbits, Saturn's gas accretion starts. 
At t X 8000 orbits, Saturn's growth depletes its coorbital re- 
gion and the positive corotation torque exerted on Saturn disap- 
pears, causing Saturn to once again migrate inward. As Saturn 
catches up with Jupiter, it is captured in the 2:1 resonance from 
t ~ 9000 to f ~ 1.3 X 10"* orbits. The resonant interaction 
causes significant growth of the planets' eccentricities up to 
ej ~ 0.3 and es ~ 0.15 (second panel of Fig. [T)- In their 
high-eccentricity state, the planets briefly repel each other once 
again outside the 2:1 resonance. Next, their eccentricities are 
quickly damped by the disk and Saturn's begins a phase of run- 
away inward migration (Masset & Papaloizou 2003). Its rapid 
migration allows Saturn to cross over the 2:1 resonance and it 
is captured in the 3:2 resonance at f ~ 1.4 x 10"* orbits. 

Once trapped in 3:2 resonance, Jupiter and Saturn migrate 
outward for a short time but, in contrast with models II and 

13, the outward migration is not maintained. Instead, the plan- 
ets remain on roughly stationary orbits with semimajor axes of 
ay ~ 1.3 and as ~ 1.7. Fig. |6] shows the evolution of the disk 
torques for the three models at times when a quasi-stationary 
state is reached. As expected, the disk torque experienced by 
Jupiter is positive in models II and 13 but it oscillates about 
zero in model 14. The total torque exerted on Saturn is clearly 
positive for both models II and 13, indicating that the posi- 
tive corotation torque due to disc material flowing from the 
outer disk across the gap is overcoming the negative differential 
Lindblad torque in these runs. However, the torque on Saturn 
is negative in model 14. 

To understand the origin of the unexpected zero torque ex- 
erted on Jupiter in model 14, Fig. Q shows the disk's surface 
density profile at f ~ 1.5 x 10"* orbits for the three models. 
Compared with runs II and 13, the surface density at the posi- 
tion of Jupiter is much higher in model 14. Indeed, the amount 
of disk material enclosed in the Hill sphere of the planets is 
about ~ 1 Me in models II and 13 and ~ 10 M® for model 

14. Fig. [8] shows that most of the mass in Jupiter's Hill sphere 
in model 13 was acquired in a short time at f « 1.2 x 10"* 
orbits, when Jupiter and Saturn's eccentricities were at their 
peak and the planets underwent large radial excursions beyond 
the edges of the gap into the gaseous disk. Although we ex- 
clude gas material located inside the planet's Hill sphere from 
the torque calculation, examination of the torque distribution 
shows that the large amount of gas material located in the vicin- 
ity of Jupiter's Hill sphere can indeed contribute significantly 
to the total torque exerted on that planet (e.g. Crida et al. 2009). 
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Fig. 8. Evolution of the disk mass located inside the Hill sphere 
of the planets for models II (black), 13 (red) and 14 (blue). 



To demonstrate that this disk region is responsible for the stop- 
ping of migration, we restarted the run 14 at f ~ 1.5 x 10"* 
orbits while slowly removing the gas within each planets' Hill 
sphere but while keeping the planets' actual masses fixed. The 
results of this calculation are shown in Fig. |9] and clearly indi- 
cate that outward migration is recovered by removing the resid- 
ual gas material from the vicinity of Jupiter Given that most of 
the gas bound to Jupiter was acquired at relatively late time, 
the "correct" outcome of model 14 still includes a complicated 
orbital evolution (including temporary capture in the 2:1 res- 
onance and eccentricity excitation) but an outward, sustained 
migration at later times once the planets are trapped in the 3:2 
resonance. The extra gas within Jupiter's Hill sphere, whose 
inertia was preventing outward migration, should realistically 
have been accreted or repelled by Jupiter on a relatively short 
time scale and should not inhibit outward migration. 

Nevertheless, we comment that the high density region lo- 
cated in the vicinity of Jupiter results from the use of an isother- 
mal equation of state, for which pressure gradients corre- 
spond to density gradients alone. In non-isothermal disk mod- 
els, this density peak would be reduced since pressure gradi- 
ents are partly supported by temperature gradients in that case 
(Paardekooper & Mellema 2008). Moreover, in these cases, it 
is not expected that such a large amount of disk material would 
accumulate in the planet's vicinity. Indeed, this would corre- 
spond to a large accretion rate onto Jupiter which would result 
to a significant heating of the Roche lobe, preventing thereby 
further accretion (Peplinski et al. 2008). 

3.2.2. Effect of the disk viscosity: model 15 

Fig.fTOlillustrates the evolution of a simulation in which a was 
set to a = 4 X lO"-' (run 15). As a consequence, both the growth 
and type II migration timescales are shorter Indeed, the accre- 
tion rates onto the planets are clearly enhanced in run 15 com- 
pared with run II and Jupiter's growth timescale is reduced by 
~ 30%. Because the viscosity is higher, we expect gap opening 
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Fig. 9. Evolution of the semimajor axes of Jupiter (solid line) 
and Saturn (dashed Une) for model 14 (black) and for a restart 
run in which gas is continuously removed from the Hill spheres 
of the planets (red). 

to occur later in the growth of Jupiter, i.e., for a higher value of 
nij. Using the criterion for gap opening derived by Crida et al. 
(2006): 

^ H 50 ^ 

where q = hij/Mq, Rh - aj(mj/3MQy^^ is the Hill radius of 
Jupiter and K - a^Q.j/v is the Reynolds number, we indeed 
predict that gap opening should occur for nij > 0.23 Mj in run 
II and for nij > 0.35 Mj in model 15. 

Thus, for a higher viscosity, Jupiter's gap grows later and 
its type II migration is faster. This means that, when Saturn's 
gas accretion starts, the Jupiter-Saturn separation is larger for 
the case of a higher viscosity. Indeed, for run 15 the two planets 
are significantly farther apart than for run II, just interior to the 
2: 1 resonance. For even higher viscosities (higher values of a), 
Saturn's core would be pushed beyond the 2:1 resonance with 
Jupiter such that subsequent evolution could involve temporary 
capture in this resonance (as in run 14 discussed above). 

In model 15, the early stages of Saturn's growth involve 
convergent migration of the two planets followed by trapping 
in the 3 : 2 resonance. Here, reversal of migration occurs for 
slightly higher planet masses than seen previously - nij = 0.8 
Mj and m, - 0.3 Mj - but the final outcome is the same, 
namely sustained outward migration with the planets maintain- 
ing their 3 : 2 commensurability. We also find that the outward 
migration is slower for higher a because Jupiter's gap becomes 
shallower as the viscosity increases, in agreement with the re- 
sults of Morbidelli & Crida (2007). 

3.2.3. Effect of the disk's aspect ratio h: models 16, 17 

We tested the effect of the disk's aspect ratio h = H/r from 
h = 0.03 (run 16) to /i = 0.05 (run 17). Fig. [lU shows the 
evolution these runs as compared with our fiducial case. From 
Eq. 3 we know that lower-mass planets can open gaps in thin- 
ner disks. Thus, for h - 0.03, Saturn's core is trapped at the 
edge of Jupiter's gap early in the simulation, and Saturn's core 
is pushed outward as Jupiter's mass increases and as its gap 
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Fig. 10. Upper (first) panel: Evolution of the planet masses for 
models II (black) and 15 (red). Second panel: Evolution of 
the semimajor axes. Third panel: Evolution of Jupiter's eccen- 
tricity. Fourth panel: Evolution of Saturn's eccentricity. Fifth 
panel: Evolution of the period ratio. 



widens; this episode of outward migration of Saturn's core is 
apparent between 1000 and 2000 orbits in Fig.\n\(h = 0.03, 
second panel). The planets' orbital separation reaches a peak 
value just outside the 2:1 resonance (see bottom panel of Fig. 
fTTT l. For this run (h - 0.03) accretion onto Saturn's core is 
switched on at f ~ 1900 orbits such that during the early stages 
of its growth, Saturn still follows Jupiter's gap through the ac- 
tion of the corotation torque. At later times, the interaction with 
the disk becomes non-linear and Saturn passes through both 
Jupiter's gap, is captured in the 3:2 resonance with Jupiter, 
and the two planets migrate outward. Because Jupiter's gap is 
deeper than in the fiducial case II, the outward migration is 
faster for model 16 (see also Morbidelli & Crida 2007). 

In run 17 the disk was thicker [h - 0.05) and this resulted in 
a different mode of evolution. As before, Saturn's core was cap- 
tured at the edge of Jupiter's gap and pushed outward as the gap 
widened, this time beyond the 2: 1 resonance with Jupiter. Once 
Saturn accreted enough gas to cancel the effect of the corotation 
torque, it became trapped in the 2:1 resonance. This is because 
Jupiter's gap is shallower for the thicker disk, causing slower 
convergent migration of the two planets. Thus, Saturn is unable 
to cross the 2: 1 resonance (see also Rein et al. 2010). Of course, 
disruption of the 2:1 resonance followed by capture in the 3:2 
resonance on longer timescales can not be ruled out. Indeed, 
Pierens & Nelson (2008) showed that, for a scenario close to 
the setup of model 17, the system is temporarily locked in the 
2:1 resonance but the resonance is broken and the planets are 
evenrually trapped in 3:2 resonance. 
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Fig. 11. Upper (first) panel: the evolution of the planet masses 
for models II (black), 16 (red) and 17 (blue). Second panel: evo- 
lution of the semimajoraxes. Thirdpanel: evolution of Jupiter's 
eccentricity. Fourth panel: evolution of Saturn's eccentricity. 
Fifth panel: evolution of the period ratio. 



3.2.4. Effect of the disk surface density profile: models 
18, 19 

We now test the effect of the disk's radial surface density pro- 
file, where the surface density S varies with orbital radius R 
as S oc R°'. We compare two runs with cr - 312 (models 18 
and 19) with our standard models that have cr - 1/2 (II and 
14). We note that a cr = 1/2 profile corresponds to disks with 
constant accretion rates and/? = 1, whereas sub-mm measure- 
ments of young protoplanetary disks appear to favor cr « 0.5 - 1 
(e.g., Mundy et al. 2000; Andrews & WilUams 2007) and dif- 
ferent interpretations of the minimum-mass solar nebula model 
(Weidenschilling 1977; Hayashi 1981) yields values of cr be- 
tween 1/2 (Davis 2005) and 2 (Desch 2007). 

Eq. [Tlpredicts that for cr - 1.5, Jupiter and Saturn's cores 
should migrate at the same rate (in an isothermal disk) rather 
than undergoing convergent migration. The consequence is 
that, in contrast with models in which cr - 1/2, the cores are 
not locked in 3:2 resonance when Jupiter starts to accrete gas 
from the disk. Fig. [12] shows the evolution of run 18; for this 
case accretion onto Saturn's core starts when Jupiter has grown 
to half of its full mass. In run 19 gas accretion is slower than 
for model II simply because the disk mass included within the 
feeding zones of the planets is smaller for larger values of cr. 

The evolution of run 18 is virtually identical to run II but 
slower Both the accretion rates and the migration rates are 
slower for cr = 3/2 (19). This is simply because of the smaller 
outer disk mass in a disk with a steep surface density profile; 
the annular mass scales as /?' °" and the mass within a planet's 
Hill sphere scales as 7?^ °". Thus, in terms of the outward migra- 



tion of Jupiter and Saturn, the mass flux across the gap is sig- 
nificantly smaller for the disk with cr - 3/2. To illustrate this. 
Fig. [T3]shows the disk surface density for runs II and 18 when 
outward migration is about to be triggered and at f = 10"* or- 
bits. The smaller gas flux across the gap for cr - 312 decreases 
the magnitude of the (positive) corotation torque as well as the 
density in the inner disk, which has the effect of weakening the 
(also positive) inner Lindblad torque exerted on Jupiter This 
is shown in the time evolution of the disk profiles (Fig. [T3T l: 
the inner disk surface density increases with time in run II due 
to the gas flowing through the gap whereas such an effect is 
marginal in model 18. 

The disk's surface density profile had only a small impact 
on the evolution of simulations II and 18, in which accretion 
onto Saturn's core started when Jupiter reached half of its fi- 
nal mass. This is also the case for runs 14 and 19, in which 
gas accretion onto Saturn's core is switched on once Jupiter is 
fully formed (Fig.[T4ll. In both runs and as discussed in Section 
3.2.1, Jupiter's eccentricity reaches values as high as ej ~ 0.3 
during capture in the 2:1 resonance. And as above, the evolu- 
tion of the simulation with a steeper disk density profile (run 
19) is slower than for the shallower profile (13). Thus, although 
the early evolution of run 19 resembles a stretched-out version 
of run 13, the final fate of the system is still not reached despite 
the very long timescale covered by the simulation (~ 3 x 10'* 
orbits). However, the fifth panel of Fig. [14]- which displays the 
time evolution of the period ratio - suggests that the 2:1 reso- 
nance will be disrupted in the next ~ lO'* orbits. Extrapolating 
the results from model II, subsequent evolution should involve 
capture in 3:2 resonance followed by outward migration of the 
Jupiter-Saturn system. 

3.2.5. Effects of the gas disk's dispersion 

As the Solar Nebula dispersed, the giant planets' migration and 
accretion stopped. We address the effect of the dispersion of 
gas disk dispersion on the simulation presented in Sect. 13.11 
using four additional simulations in which, after a delay, the 
gas surface density was forced to decay exponentially with an 
e-folding time Tdisp- We used model parameters as in run II 
and we assumed that both the gas disk dispersion and accretion 
onto Saturn's core start at the same time. In these simulations, 
we varied the value of r^z/jp which was set to r^usp = 10-^, 3 x 10^ 
and 10"* orbits respectively. 

Fig. [15] shows the evolution of simulations with different 
values of T^isp. The run with Tjisp = 10-' orbits is clearly not a 
viable scenario since the gas lifetime is so short that accretion 
is cut off early and Jupiter and Saturn never reach their true 
masses. For Tdtsp = 3 x 10"^ and T^isp = 10"* orbits, however, 
the correct masses for Jupiter and Saturn are obtained and the 
system migrates outward, reaching final orbits of aj ~ 1 .8 and 
as ~ 2.4 for Tdisp = 3 x 10^ and aj ~ 3.5 and as ~ 6.5 for 
Tdisp = 10"* orbits. 

The series of simulations presented in this section suggest 
that a higher value for Tdisp is required for such a model to be 
consistent with the "Grand Tack" scenario. An alternate possi- 
bility is that both Jupiter and Saturn formed early in the lifetime 
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Fig. 12. Upper (first) panel: the evolution of the planet masses 
for models II (black) and 18 (red). Second panel: evolution of 
the semimajor axes. Third panel: evolution of Jupiter's eccen- 
tricity. Fourth panel: evolution of Saturn's eccentricity. Fifth 
panel: evolution of the period ratio. 





Fig. 14. Upper (first) panel: the evolution of the planet masses 
for models 11 (black) and 19 (red). Second panel: evolution of 
the semimajor axes. Third panel: evolution of Jupiter's eccen- 
tricity. Fourth panel: evolution of Saturn's eccentricity. Fifth 
panel: evolution of the period ratio. 
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Fig. 13. Upper panel: the disk surface density profile for model 
11 prior that outward migration of the Jupiter and Saturn sys- 
tem occurs (solid line) and at the end of the simulation (dashed 
line). Lower panel: same but for model 18 



of the Solar Nebula, long before the disk was being dispersed. 
To investigate this question, we performed an additional run 
in which the gas disk disperses when Jupiter and Saturn ap- 
proach their current orbits. Fig. [16] shows the results of a sim- 
ulation with Tdisp - 1000 orbits and in which disk dispersion 
was initiated after tdisp ~ 2.2 x 10"^ orbits. As expected, the ec- 
centricities grow due to the disk induced eccentricity damping 
being cancelled and saturate at ej ~ 0.05 and es ~ 0.1. This 
effect is also compounded by the fact that the planets become 
locked deeper in the 3:2 resonance while the gas is being dis- 
persed. Here, the planets reach final orbits with aj ~ 3.6 and 
as ~ 4.7 but it is clear that a similar simulation performed with 
an adequate value for tacc would lead to both Jupiter and Saturn 
reaching their expected pre-Nice model orbits with aj * 5.4 
AU (Tsiganis et al. 2005). 

4. Simulations in radiative disks 

To estimate the influence of a more realistic treatment of the 
disk thermodynamics, we performed two additional simula- 
tions that do not use a locally isothermal equation of state but 
in which the full energy equation is solved (labelled as Rl and 
R2). For run Rl, the parameters for the disk and planets were 
the same as in model 11 whereas for run R2, parameters were 
identical to run 13. 

For these radiative disk simulations, we ran a preliminary 
model without any planet in order to obtain a new equilibrium 
state for the disk where viscous heating balances radiative cool- 
ing from the disk surfaces. We then restarted the simulation 
with the cores of Jupiter and Saturn embedded in the disk. Fig. 



Pierens & Raymond: Two-phase migration of Jupiter and Saturn 



11 



1 0000 orbits 




2000 4000 6000 8000 10000 
Time (orbits) 

Fig. 15. Upper (first) panel: the evolution of the planet masses 
for simulations in which gas disk dispersion is considered. 
Second panel: evolution of the semimajor axes. Third panel: 
evolution of Jupiter's eccentricity. Fourth panel: evolution of 
Saturn's eccentricity. Fifth panel: evolution of the period ratio. 
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Fig. 16. Upper left (first) panel: evolution of the planet masses 
for a simulation with t^, 



sp 



1000 orbits and in which onset 
of disk dispersion occurs at f ~ 2.2 x 10"^ orbits. Upper right 
(second) panel: evolution of the semimajor axes. Third panel: 
evolution of the eccentricities. Fourth panel: evolution of the 
period ratio. 



[TTl shows the density and temperature profiles when a station- 
ary state for the disk is reached. Although the surface density at 
the initial positions of Jupiter (aj - 2) and Saturn (as - 2.65) 
is similar between the isothermal and radiative disks, the tem- 
perature at these locations are somewhat lower in the radiative 
calculation. Thus, the disk's aspect ratio is HIR ~ 3.3 x 10"^ 
at the initial position of Jupiter and H/R ~ 2.8 X 10"^ at the 
initial location of Saturn. For the disk model considered here, 
we note that the disk is initially optically thick for R < 6. 

Fig. [T8] compares the evolution of the system for both the 
isothermal model II and the radiative calculation Rl . Although 



including heating/cooling effects appears to have little impact 
on the mass-growth history of Jupiter and Saturn, the dynami- 
cal evolution is very different between the isothermal and radia- 
tive disks. For model Rl, the 10 M® cores of Jupiter and Saturn 
initially migrate much more slowly than for run II . This occurs 
because the entropy gradient within the horseshoe region of 
the planets gives rise to a positive corotation torque (Baruteau 
& Masset 2008; Paardekooper & Papaloizou 2008) that acts 
in opposition to the negative differential Lindbald torque. This 
positive corotation torque can be sustained provided that diffu- 
sive processes (thermal diffusion, heating/cooling effects...) can 
restore the original temperature profile and that the diffusion 
timescale across the horseshoe region is shorter than the libra- 
tion timescale Tut- In the simulations presented here, the dif- 
fusion timescale corresponds to the vertical cooling timescale 
Tcooi = CylLT/Q where Q is the local radiative cooling and 
c,, the specific heat at constant volume. For Jupiter, Tm, ~ 42 
Torb where Tgrb is orbital period of the planet and TcooI ~ 280 
Torb while for Saturn Tub ~ 56 Torh and TcooI ~ 251 Torb, 
which means that the corotation torque is partially saturated 
for both embryos. As Jupiter grows, this corotation torque be- 
comes strong enough to push the planet outward, which is ap- 
parent at f ~ 1000 orbits in the second panel of Fig. [18] When 
Jupiter's mass has reached ~ 45 M® however, the planet opens 
a gap in the disk which consequently suppresses the corotation 
torque and makes Jupiter migrate inward again on the type II 
migration timescale. 

Jupiter's type II migration is faster than Saturn's type I mi- 
gration so the two planets' orbits diverge. This can be seen for 
t < 2000 orbits in the fifth panel of Fig. [18] in this case, the 
period ratio increased to {aslajY'^ ~ 2.4 before Saturn opened 
a gap in the disk. From this point in time, Jupiter and Saturn 
migrated convergently until they became captured in the 2:1 
resonance. The planets then migrated outward for a brief in- 
terval but the 2: 1 resonance configuration became unstable and 
was broken at f ~ 1.2 x 10"^ orbits. Jupiter and Saturn then be- 
came temporarily trapped in 5:3 resonance, during which time 
slow outward migration continued, but once again the reso- 
nance was broken. Next, Saturn became locked in the 3:2 res- 
onance with Jupiter and the two planets migrated outward to- 
gether. However, in contrast with the isothermal runs, this con- 
figuration proved unstable, as Jupiter and Saturn underwent a 
dynamical instability leading to a weak scattering event that 
launched Saturn beyond the 2: 1 resonance at f ~ 2 x 10"* orbits. 
At the end of the run, the final fate of the run is still uncertain 
- it is possible that the subsequent evolution will again involve 
a cycle of temporary capture in the 2:1, 5:3 and 3:2 resonances 
followed by instabilities until the disk dissipates. 

We think that this instability arose because in the outer parts 
of the radiative disk both the temperature and the disk aspect 
ratio decrease. Such a cold disk acts both to accelerate outward 
migration and to decrease disk-induced eccentricity damping. 
Thus, as the two planets migrated outward their eccentrici- 
ties were significantly higher than for the isothermal runs (see 
the third and fourth panels of Fig. \T% . It is well-known that, 
for the case of two inward-migrating planets, the inner disk 
mass plays a key role in damping the inner planet's eccentric- 
ity and thus maintaining dynamical stability (Crida et al. 2008). 
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Because the disk aspect ratio is smaller, the gap edge here lies 
further from the planets compared with isothermal runs, re- 
sulting in a weaker eccentricty damping from the inner disk. 
Simulation Rl appears to present a similar scenario but with 
two outward-migrating planets, with Saturn's relatively large 
and chaotically-varying eccentricity acting as the trigger for 
instability. It should be noted that the outer parts of the ra- 
diative disk are so cold because viscosity is the only heating 
process in the simulation. However, it is well known that stel- 
lar irradiation is the main heat source at R > 3 AU (D'alessio 
et al. 1998). Thus, a more realistic disk should probably have 
a warmer outer disk. It is unclear if this would encourage out- 
ward migration by reducing the likelihood of instability or dis- 
courage outward migration by overly puffing up the disk. This 
is an area for future study. 

Fig. [19] show the evolution of radiative run R2 in which 
Jupiter and Saturn start to accrete gas at the same time. As it 
grows, Jupiter's outward migration due to the entropy-related 
corotation torque makes the planets converge, until their orbital 
period ratio reaches a minimum of ~ 1 .25 at f ~ 800 orbits (i.e., 
Saturn is interior to the 3:2 resonance). Once it has accreted 
enough gas to cancel the effect of the corotation torque, Jupiter 
migrates inward again, resulting in a divergent migration which 
continues until the planets become locked in the 3:2 resonance. 
Outward migration of Jupiter and Saturn is then triggered and 
appears to be maintained until the simulation was stopped after 
10^ orbits. Compared with the isothermal disk, Saturn's eccen- 
tricity is significantly higher although its behavior is steady and 
not obviously chaotic as in run Rl . We do not know if this out- 
ward migration will continue indefinitely or whether the system 
might be subject to an instability similar to run Rl. 

Thus, simulations Rl and R2 demonstrate that periods of 
outward migration of Jupiter and Saturn in radiative disks are 
viable. However, only one of two simulations produced a clear 
two-phase migration. Given the limitations in our simulations 
(especially with regards to the thermal state of the outer disk), 
we do not know whether a two-phase migration of Jupiter and 
Saturn is a likely outcome in radiative disks. Indeed, these ra- 
diative simulations were performed assuming that the radius 
/? = 1 in the computational domain corresponds to 5 AU. 
Contrary to isothermal runs, it is worth to note that results 
from radiative simulations can not be scaled to apply to dif- 
ferent parameters. Unless the Grand Tack occured in the last 
stages of the disk's lifetime, our radiative calculations there- 
fore probably underestimate the disk temperature at the lo- 
cation where Jupiter's migration reversed. We are currently 
working to test the outward migration mechanism of Masset 
& Snellgrove (2001) in more realistic radiative disks under a 
range of physical conditions. 



5. Evolution of the Solar Nebula 

Our results thus far show that a two-phase migration of Jupiter 
and Saturn is extremely robust in isothermal disks but is as-yet 
uncertain in radiative disks. Protoplanetary disks can be con- 
sidered to be isothermal if they are optically thin (i.e., if the 
optical depth r < 1) and radiative if they are optically thick 




Isothermal 
.Radiative 



Fig. 17. Surface density (upper panel) and temperature (lower 
panel) profiles at equilibrium for the isothermal (black) and ra- 
diative (red) models. 

(r > 1). But when in the Solar Nebula's history was it isother- 
mal or radiative? 

To address this question we constructed a simple, 1-D 
model of the viscously-evolving Solar Nebula. The disk ex- 
tended from 0.1 to 40 AU and initially contained 40 Mj fol- 
lowing an /? '^^ surface density profile. We adopted an a pre- 
scription for the disk's viscosity (Shakura & Sunyaev 1973) 
and used the same value as in most of the hydro simulations, 
a - 2 X lO"-'. We solved the viscous diffusion for the surface 
density: 



'dt 



3_d_ 
RdR 



^/R 



dy/RvY 



dR 



(4) 



and calculated the temperature using a simple radiative balance 
between the disk's viscous heating and radiative cooling (as in 
Lyra et al 2010): 



lo-T'^ ^T,ff\-vl.Q} 



(5) 



where <j is the Stephan-Boltzmann constant and Q is the orbital 
frequency. The effective optical depth to the disk midplane is 
represented by t^//, which is defined as 



3t 



V3 l_ 

4 ^ At' 



(6) 



The optical depth is t = /cS/2. We assume that the opacity k is 
dominated by small grains and use the values from Bell & Lin 
(1994). 

Fig. |20] shows 8 Myr in the evolution of a representative 
Solar Nebula. As the disk viscously spreads its surface density 
decreases uniformly and the disk cools. The cooling is not uni- 
form due to the large variations in opacity between different 
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Fig. 18. Upper (first) panel: the evolution of the planet masses 
for the isothermal model II (black) and the radiative model Rl 
(red). Second panel: evolution of the semimajor axes. Third 
panel: evolution of Jupiter's eccentricity. Fourth panel: evolu- 
tion of Saturn's eccentricity. Fifth panel: evolution of the period 
ratio. 
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Fig. 19. Upper (first) panel: the evolution of the planet masses 
for the isothermal model 13 (black) and the radiative model R2 
(red). Second panel: evolution of the semimajor axes. Third 
panel: evolution of Jupiter's eccentricity. Fourth panel: evolu- 
tion of Saturn's eccentricity. Fifth panel: evolution of the period 
ratio. 



temperature regimes. Similarly, the disk is initially optically 
thick in its inner 20 AU and optically thin farther out. In time, 
the boundary between optically thick and thin moves inward 
but interior to 1 AU the disk remains optically thick throughout. 
Because we have not included photo-evaporation, the disk's 
density continues to decrease but never to zero. In reality, at 
some point we expect the disk to be completely removed by 
either photo-evaporation (Hollenbach et al 1994; Adams et al. 
2004) or perhaps an MRI-related instability (Chiang & Murray- 
Clay 2007). The optical depth in a given location depends on 
the local disk properties, although the evolution of the disk's 
surface density profile certainly depends on whether the disk 
dissipates from the outside-in or the inside-out. 

How does the Grand Tack fit in the context of this simple 
model? The region of interest, from roughly 1-10 AU, is clearly 
in the radiative regime early in the disk's lifetime. In the last 
few Myr this region transitions to an isothermal state. The rel- 
evant boundary where an isothermal disk should affect Jupiter 
and Saturn's evolution is the outer edge of Saturn's gap when 
Saturn is at its closest to the Sun. This corresponds to about 2.5 
AU for Jupiter at 1.5 AU and Saturn in 2:3 MMR at 1.97 AU. 
At 2.5 AU, the disk transitions from radiative to isothermal af- 
ter roughly 6 Myr of evolution, although the time at which this 
state is reached is parameter-dependent. 

When the disk transitions to an isothermal state its density 
is significantly decreased compared to its initial state or to the 
configuration of the hydrodynamical simulations presented in 
Sects. 3 and 4. Could Jupiter and Saturn migrate from 1 .5-2 AU 
out to their current locations in such a low-mass disk? There 



are two criteria that should be required to allow for long-range 
outward migration via the Masset & Snellgrove (2001) mecha- 
nism. First, the inner lindblad torque acting on Jupiter must be 
larger than the outer lindblad torque acting on Saturn. Second, 
the angular momentum content of the gas though which Jupiter 
and Saturn will migrate must be sufficient to transport them a 
long distance. At the time of the radiative-to-isothermal transi- 
tion the torque balance criterion is met if we assume that the 
disk profile at that time should have been sculpted by Jupiter 
and Saturn as in the hydrodynamical simulations. In addition, 
the angular momentum content in the gas from 2.6 to roughly 
8-10 AU is ~ 1.5 times larger than that needed to move Jupiter 
and Saturn to 5.4 and 7.1 AU. 

Despite meeting the theoretical criteria for outward migra- 
tion, we ran two additional hydrodynamical simulations to see 
if outward migration of Jupiter and Saturn could truly occur in 
such low-mass disks. One simulation was run with an isother- 
mal equation of state whereas the other included radiative ef- 
fects. Jupiter and Saturn started the simulations fully-formed 
and were placed just exterior to the 3:2 resonance. The disk 
mass interior to Jupiter's orbit was only ~ 0.4My.In addition, 
in these two simulations the spatial units were AU rather than 
multiples of 5 AU (and the corresponding time units years 
rather than 5-^^^ ~ 10 years) such that these simulations truly 
test the Grand Tack at its correct scale. Given the large compu- 
tational expense these simulations were only run for 1000 years 
and serve mainly as a proof of concept. 

Fig. |2T] shows the time evolution of the semimajor axes 
of Jupiter and Saturn in these two simulations. As predicted, 
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Fig. 20. Upper (first) panel: Surface density profiles for the 
one-dimensional disk model from t - \ Myr (upper curve) to 
/ = 8 Myr. Second panel: Evolution of the temperature. Third 
panel: Evolution of the optical depth. 



Saturn became trapped in resonance and the two planets tacked 
and migrated outward in both cases. This confirms the results 
found in previous sections and shows that the Grand Tack 
mechanism applies on the relevant spatial scale. 

It therefore appears that outward migration of Jupiter and 
Saturn from 1.5-2 AU to beyond 5 AU is a natural outcome in 
an isothermal Solar Nebula. Of course, in certain situations out- 
ward migration may occur in a radiative disk (e.g., FigfTSl). But, 
for the limiting case in which outward migration can only occur 
in an isothermal disk, we still expect the Grand Tack to happen 
because in the last stages of the disk's lifetime it is necessarily 
optically thin and, as seen in Sect. 3, this leads inevitably to 
outward migration. 

6. Summary and Discussion 

Our results indicate that Jupiter and Saturn probably under- 
went a two-phase, inward-then-outward migration. In our sim- 
ulations, Jupiter and Saturn start as 10 M^ cores and type I 
migrate; inward for isothermal disks, inward or outward for ra- 
diative disks. In most cases the two cores become locked in 
3:2 mean motion resonance (MMR). At this point or after a 
delay of 500-1000 orbits, we allowed Jupiter to start accret- 
ing gas from the disk. When Jupiter reaches the gap-opening 
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Fig. 21. Time evolution of the semimajor axes of Jupiter (lower 
line) and Saturn (upper line) for hydrodynamical simulations 
adapted to the Solar Nebula. 



mass, it undergoes a phase of rapid inward migration as it 
clears out its gap (sometimes called type III migration; Masset 
& Papaloizou 2003) then settles into standard, type 11 migra- 
tion. Inward migration continues until Saturn accretes enough 
gas to reach the gap-opening mass itself. At this point, Saturn's 
inward migration accelerates and is again trapped in the 3:2 
MMR with Jupiter Outward migration of both giant planets 
is then triggered via the mechanism of Masset & Snellgrove 
(2001). Outward migration stops when either a) the disk dissi- 
pates (as in Sect. 13.2.5b . b) Saturn reaches the outer edge of the 
disk, or, c) if the disk is flared, the giant planets drop below the 
local gap-opening mass (e.g., Crida et al. 2009). 

An additional stopping - or at least slowing - mechanism 
exists if the planets are unable to maintain a well-aligned res- 
onant lock during migration. For example, in simulation II 
Jupiter and Saturn's rate of outward migration slowed signif- 
icantly when one resonant angle transitioned from libration to 
circulation (see Figs. 1 and 3). Here, this appears to be due to 
a positive corotation torque exerted on Saturn by gas that pol- 
luted Jupiter and Saturn's common gap as the distance between 
the two planets increased during outward migration. On longer 
timescales it is unclear if this mechanism would continue to 
slow down and eventually stop the outward migration. 

In isothermal disks, the two phase migration of Jupiter 
and Saturn holds for almost the full range of parameters that 
we tested (Sect. 3). The only situation for which this result 
does not hold is if the gaseous Solar Nebula is relatively thick 
{h - H/r > 0.05). Both the disk's surface density profile and 
the value for the disk aspect ratio had an effect on the migration 
rate: disks with either shallower profiles or lower values of h re- 
sult in faster migration. Changing the disk's viscosity had little 
effect on the outcome, although we only tested a very small 
range. In one simulation (14; Sect. 13.2.11 ). Saturn's core was 
pushed past the 2:1 MMR with Jupiter leading to a significant 
eccentricity increase for both planets before the resonance was 
crossed, Saturn was trapped in the 2:3 MMR and both planets 
migrated outward. This dynamic phase of resonance crossing 
and eccentricity excitation is likely to be quite sensitive to the 
detailed properties of the disk (e.g., the scale height and viscos- 
ity) that determine the gap profile. 
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We performed two simulations in radiative disks with 
mixed outcomes (Sect. 4). In the first case, Jupiter and Saturn 
started accreting together and so stayed relatively close to each 
other The planets became locked in the 3:2 MMR and migrated 
outward even faster than in isothermal simulations due to the 
small aspect ratio of the outer disk (h x 0.03). In the second 
case, Saturn started to accrete when Jupiter reached half its fi- 
nal mass (i.e., xj = 0.5), by which time the two planets were 
beyond the 2: 1 MMR. The planets succeeded in breaking the 
2:1 and 5:3 MMRs, became trapped in the 3:2 MMR and mi- 
grated outward only to undergo a dynamical instability putting 
the planets once again beyond the 2:1 MMR. A more detailed 
study of outward migration in radiative disks is underway. 

Using a simple 1-D model of an evolving Solar Nebula we 
showed that the disk should be optically thick at early times, 
then transition to optically thin from the outside-in during the 
late phases of its evolution. At the orbital distance in question 
(1-10 AU), the disk transitions from radiative to isothermal be- 
havior in the last 1-2 Myr of its evolution. Thus, even if we 
make the "pessimistic" assumption that an isothermal disk is 
required for outward migration of Jupiter and Saturn, the disk 
fulfills the criteria for long-range outward migration in its late 
phases. Outward migration of Jupiter and Saturn at this time is 
very likely provided the disk remains thin (h < 0.05). 

Our simulations therefore show that an inward-then- 
outward migration of Jupiter and Saturn is extremely likely, 
and that the last phase of outward migration probably coincided 
with the late phases of the dissipation of the Solar Nebula. 
This is of particular interest because the two phase migration 
of Jupiter and Saturn helps resolve a long-standing problem in 
terrestrial planet formation. For over 20 years, simulations of 
terrestrial accretion have been unable to reproduce Mars' rel- 
atively small mass (0.11 M®; Wetherill 1978, 1991; Chambers 
2001; Raymond et al. 2009). This problem arises because, in 
a Solar Nebula that varies smoothly in orbital radius, there is 
a comparable or larger amount of mass in the vicinity of Mars 
than the Earth. For Mars to be so much smaller than Earth, 
most of the mass between roughly 1-3 AU must be removed 
(e.g., Raymond et al. 2006, 2009; O'Brien et al. 2006). Several 
mechanisms have been proposed to remove this mass, includ- 
ing strong secular resonances (Thommes et al. 2008, Raymond 
et al. 2009) and a narrow dip in the surface density caused by a 
radial dependence of the disk's viscosity (i.e., a dead zone; Jin 
et al. 2008). However, the problem is most easily and much bet- 
ter solved if the terrestrial planets did not form from a wide disk 
of planetary embryos but instead from a narrow annulus ex- 
tending only from 0.7-1 AU (Wetherill 1978; Chambers 2001; 
Hansen 2009). In that case. Mars' small mass is simply an edge 
effect: Mars is small was built from one or perhaps a few em- 
bryos that were scattered beyond the edge of the embryo disk 
(this is also the case for Mercury, which was scattered inward 
beyond the inner edge of the embryo disk). In contrast. Earth 
and Venus formed within the annulus and are consequently 
much more massive. Simulations of terrestrial planet formation 
can quantitatively reproduce the orbits and masses of all four 
terrestrial planets as well as their radial distribution (Hansen 
2009). 



The flaw in simulations of terrestrial planet formation in 
truncated disks is that they had no justification for the trunca- 
tion; the ad-hoc initial conditions were simply chosen because 
they provided a good fit to the actual terrestrial planets (Hansen 
2009). The two phase migration of Jupiter and Saturn provides 
such a justification via the Grand Tack model of Walsh et al. 
(2011). If Jupiter's turnaround point was at ~ 1.5 AU then 
it would have naturally truncated the inner disk of embryos 
and planetesimals at about 1 AU - in most of our simulations 
Jupiter indeed tacked at roughly this distance. As expected, the 
terrestrial planets that form from this disk quantitatively re- 
produce the actual terrestrial planets (Walsh et al. 2011). The 
Grand Tack model also provides the best explanation to date 
for the observed dichotomy between the inner and outer aster- 
oid belt (Gradie & Tedesco 1982). Thus, the present-day Solar 
System appears to bear the imprint of a two phase migration of 
Jupiter and Saturn. Our hydrodynamical simulations provide 
support for the Grand Tack scenario. 

As with any numerical study, our simulations do not fully 
represent reality. The aspect of our simulations that is proba- 
bly the least realistic is the gas accretion onto the giant planets' 
cores. In our simulations, gas accretion onto Jupiter and Saturn 
is extremely fast. Once accretion starts, Jupiter and Saturn 
reach their final masses in only a few thousand years, whereas 
the Kelvin-Helmholtz time in protoplanetary disks is more like 
~ 10^ years. In addition, accretion onto growing giant planet 
cores requires transferring gas through circum-planetary ac- 
cretion disks whose physical properties are poorly constrained 
(e.g.. Ward & Canup 2010). Once the planets reached their ac- 
tual masses we artificially turned off gas accretion. If, during 
the outward migration Saturn accreted enough gas to carve a 
gap as deep as Jupiter's then Saturn's outer lindblad torque 
would balance Jupiter's inner lindblad torque, outward migra- 
tion would stop and the planets would turn back around and mi- 
grate inward. The impact of a more realistic accretion history 
on Jupiter and Saturn's migration remains an open question, in 
particular with regards to the interplay between gas accretion 
and the dispersal of both the cimcumstellar and cicumplane- 
tary disks. 
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